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^ , Abstract. Entropy accumulation near a quantum critical point was expected based on general 

S_i ■ scaling arguments, and has recently been explicitly observed. We explore this issue further in 

' two canonical models for quantum criticality, with particular attention paid to the potential 

effects beyond hyperscaling. In the case of a one-dimensional transverse field Ising model, we 
' derive the specific scaling form of the free energy. It follows from this scaling form that the 

singular temperature dependence at the critical field has a vanishing prefactor but the singular 
field dependence at zero temperature is realized. For the spin-density-wave model above its 
' O ' upper critical dimension, we show that the dangerously irrelevant quartic coupling comes into 

the free energy in a delicate way but in the end yields only subleading contributions beyond 
hyperscaling. We conclude that entropy accumulation near quantum critical point is a robust 
property of both models. 



1. Introduction 
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t^J- ■ Quantum critical points (QCPs) have been extensively studied in heavy fermion metals and 

related systems [I]. They occur as a non-thermal control parameter is tuned to a second-order 
phase transition at zero temperature. For thermodynamics, it was shown [2] based on scaling 
considerations that the thermal expansion [a = (l/V)(dV / dT) p ^ oc dS/dp, the variation of 
entropy S with pressure p] is more singular than the specific heat [c p = (T/N)(dS/dT) p ], in 
the general case when the tuning parameter is linearly coupled to pressure. Correspondingly, 
the Griineisen ratio, T = a/c p , diverges at the QCP, and the entropy is maximized. The 
same conclusions apply to a field-tuned QCP, where the magnetic Griineisen ratio is the 
magnetocaloric effect. The predicted divergence of the Griineisen ratio is by now widely observed 
in quantum critical heavy- fermion metals [3j 0] , and the entropy enhancement has recently been 
explicitly observed in a quantum critical ruthenate [5]. 

The scaling arguments proceed as follows. Near a QCP, the critical part of the free energy 
takes the hyperscaling form F = FoT d/ ' z+1 f '(r /T 1 ^ uz " > ) , where r = p — p c , d is the spatial 
dimension and u, z are respectively the correlation-length and dynamic exponents. The universal 
function f{x) has different asymptotic behaviors in the x — > and x — > ±oo limits, respectively 
corresponding to the quantum critical and quantum disordered/renormalized classical regimes. 
The divergence of T can be readily derived, in universal forms as 1/r for |r| 3> T and as 1/T 1 ^ UZ 
for |r| <C T [2]. The 1/r divergence in the low-temperature limit also amounts to a sign change 
of r across the QCP. As seen in Fig. [H this corresponds to an entropy maximization at the QCP 



in the low-temperature limit. The extension of this sign change to the finite temperature phase 
transitions has also been discussed [6]. 

The hyperscaling form for the free energy is based on the existence of a single critical energy 
scale, i.e., the gap for the quantum critical excitations A ~ r uz . In this paper we explore how 
the scaling form of the free energy arises in some specific models for QCPs, paying particular 
attention to the possible effects that go beyond hyperscaling. 




Figure 1. Divergence of the Griineisen ratio and the accumulation of entropy near QCPs. 



2. One-dimensional transverse field Ising model 

Consider first the one-dimensional transverse field Ising model (ID TFIM), defined by the 
Hamiltonian [7]: 

i 

where J > is the ferromagnetic Ising coupling between the neighboring spins and gJ is the 
transverse magnetic field. With the tuning of g, there exists a quantum phase transition from a 
ferromagnet to quantum paramagnet at g c = 1 [TJ [H] . Using a Jordan- Wigner transformation, we 
can represent the Hamiltonian in the thermodynamic limit in terms of free fermions and the free 
energy can be written as F = —NksT [in 2 + ~ f^d k In cosh.^^/ 2k bT)^, with the dispersion 

1/9 

= 2 J (l + g 2 — 2g cos A;) . The normalized free energy density is 

/ = = fo(g) - l -[ dk ln I 1 + e ' 2A ) : A = V(l-2) 2 A 2 + 4(5A 2 )sin 2 (fc/2) (2) 

where t = /c^T/J, and fo(g) is the ground state energy. Our focus will be on the T-dependent 
part, fi(g,T) = f(g,T) - f (g). 

Consider first the quantum critical regime, where \g — 1| /t <C 1. At temperatures small 
compared to J, we can introduce a momentum cutoff k c ~ t, below which |sin (k/2)\ /t 1 and 
therefore A <C 1. The free energy can be approximated by the contributions from k < k c , which 



is 
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While Eq. ([3]) contains approximations for the regular pieces, it reveals an important point. 
Compared with the general hyperscaling form, this expression is particular in that the linear 
term |1 — g\/t in the expansion vanishes. As a result, the magnetic Griineisen ratio as a function 
of temperature at the QCP stays finite. 

Consider next the low-temperature regions where |1 — g\ jt S> 1. Here A > 1 and we can 
expand ln(l + e~ 2A ) as e~ 2A and obtain 
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We can further replace by \j ^fg~ c without missing any singularity. The free energy has the 

exact hyperscaling form with z = 1. We obtain the magnetic Griineisen ratio to be 

(5) 
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which is indeed divergent at the QCP, g c = 1. The accumulation of entropy immediately follows 

(Fig. up. 

3. Spin-density-wave quantum critical point 

Quantum phase transitions in itinerant magnets are traditionally described in terms of a T = 
spin-density-wave (SDW) transition and formulated as a quantum Ginzburg-Landau theory 
with dynamic exponent z > 1 |9j. For concreteness, we take the 3D antiferromagnetic case as 
an example, i.e., with d = 3 and z = 2. The effective dimension of the Ginzburg-Landau theory 
is d + z, which is above the upper critical dimension 4. Correspondingly, the quartic coupling u, 
appearing in the action as u f 4 , is irrelevant in the renormalization-group (RG) sense. Under 
the RG transformation, u renormalizes to zero as the fixed point is reached, leaving only the 
Gaussian (quadratic) part. However, u is dangerously irrelevant [10l[TT]. For instance, in the 
quantum critical regime, the correlation length is determine as ~ r(T) = r + cuT 3 ^ 2 (where c 
is a constant, 4(n + 2)(l/3 — In 2\/2/(37r 3 )), with n being the number of components of the order 
parameter). One can consider u as introducing a new energy scale, and the scaling functions for 
generic physical quantities are expressed in terms of two variables f(r/T,uT 3 ^ 2 /r). 

Within the RG approach |10j . it is natural to focus on the free energy associated with the 
Gaussian term, Fq. The expression for Fq in the RG approach can be shown to be equal to 
a Gaussian free energy with r replaced by renormalized r(T). We will focus on the quantum 

critical regime, and write Fq = Fq^ + F^ , where 



^ = -»V [ A A I'* ~ fcoth 4- - l) tan- 1 j^fyj^ , (6) 
Jo (2vr) 3 J 2vr V 2T J r(T) + (q/A) 2 V ; 



and 



p(2)_ nV f A d*g [ r « de , e/T q 

F ° -~ nV k Wfk ^ tan W) + (.A) 2 • (7) 

Here, F q = Tog 2-2 , while A and Tq are respectively the ultraviolet momentum and 

(2) 

energy cutoffs. The leading temperature-dependent contribution from Fq turns out to be 
ny 4 ^ 3 3 rn + In 2 - ^ ln(V2 + 1)) c(u/T 3 /2 )T 3 / 2 , which is more singular than that from Fjp . 



This Gaussian form, however, is incomplete. The quartic coupling also introduces an 
explicitly linear in u term (t = T/Tq), which takes the form 



F. = -2„(„ + 2)NT uf ^ jf ' § coth ^ + ( J. + ^ ) /. *, . 

Iu= 2^J "7 C ° 2te 2 * e 2 + (re 2 * + l) 2 + w J J^f 2te 2x 1 + {re 2x + g 2 ) 2 ' ( ' 
We find that, to the linear order in u, the leading temperature-dependent contribution from F u 

(2) 

exactly cancels that of FX ' , leaving the singular contribution to the total free energy to be 
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which is given by 



T 2 2ttT 5 / 2 / 9r + cuT 3 / 2 /Tl /2 \ , s 

We can trace the singular terms included in Eq. ([6]) to the contributions of the low-energy 
degrees of freedom, and those included in Eq. to a unphysical origin from the degress of 
freedom at high energies and short distances near the cutoff scales. Eqs. ()6|9p are therefore 
expected to be the regularization prescription that is valid to all orders in u. The explicit result 
given in Eq. (|10p is consistent with those often quoted in the literature. The T 2 term represents 
the background Fermi liquid contribution. The critical part of the free energy is consistent with 
the scaling form T 5//2 (l + r/T). The contribution of the renormalized quartic coupling uT 3 ^ 2 to 
the free energy is subleading; the dangerously irrelevant parameter does not modify the leading 
singular part of the free energy from its hyperscaling form, thereby preserving the divergence of 
the Griineisen ratio and the accompanied entropy accumulation. 



4. Summary 

We have considered the thermodynamic properties of two models of quantum criticality. We 
have analyzed the potential effects beyond the most general hyperscaling forms in both models, 
and established that entropy accumulation is a robust property of both models. 
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